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ABSTRACT 



. Aims. Cosmic shear, the distortion of images of distant sources by the tidal gravitational field of the large-scale matter 
t^J- distribution in the Universe, is one of the most powerful cosmological probes. The measured shear field may be due not only 
QO to the gravitational lensing effect, but may contain systematic effects from the measurement process or intrinsic alignment 
(""*) ■ of galaxy shapes. One of the main probes for these systematics are the division of the shear field into E- and B-mode shear, 
where lensing only produces the former. As shown in a recent note, all currently used E-/B-mode separation methods for the 
, shear correlation functions £± require them to be measured to arbitrarily small and/or large separations which is of course not 
feasible in practice. 

Methods. We derive second-order shear statistics which provide a clean separation into E- and B-modes from measurements of 
C"| £±(#) over a finite interval only. We call these new statistics the circle and ring statistics, respectively; the latter is obtained 
O i 1 by an integral over the former. The mathematical properties of these new shear statistics are obtained, as well as specific 
expressions for applying them to observed data. 

Results. It is shown that an E-/B-mode separation can be performed on measurements of £± over a finite interval in angular 
separation, using the ring statistics. We furthermore generalize this result to derive the most general class of second-order shear 
i-^ statistics which provide a separation of E- and B-mode shear on a given angular interval $ m in < # < i) mas . In view of these 
• • ■ generalization, we discuss the aperture dispersion and the shear dispersion and their relation to the shear correlation functions. 
, Our results will be of practical use particularly for future cosmic shear surveys where highly precise measurements of the shear 
will become available and where control of systematics will be mandatory. 
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1. Introduction 

The weak gravitational lensing effect by the large-scale 
matter distribution in the Universe is recognized as one of 
the most powerful cosmological probes (see Mellier 1999: 
Bartelmann & Schneider 2001; Refregier 2003; Schneider 
2006 for reviews). First observed in 2000 (Bacon et al. 
2000; Kaiser et al. 2000; van Waerbeke et al. 2000; 
Wittman et al. 2000), this cosmic shear effect has ad- 
vanced tremendously over the past few years (see, e.g., 
van Waerbeke et al. 2001, 2002, 2005; Maoli et al. 2001; 
Jarvis et al. 2003, 2005; Hoekstra et al. 2002, 2005; Bacon 
et al. 2003; Semboloni et al. 2005), and large ongoing and 
planned surveys will turn cosmic shear into an important 
tool for cosmology. 

Send offprint requests to: P. Schneider 
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In order for satisfying its promises, systematic errors in 
cosmic shear measurements must be carefully controlled. 
Quite a few of those possible systematics have been iden- 
tified, including the uncertainties in the overall shear cal- 
ibration resulting from the correction of image smearing 
by seeing (e.g., Huterer at al. 2006) and uncertainties in 
the redshift distribution of source galaxies (e.g., Ma et al. 
2006). 

One of the important checks for possible systematic 
effects consists in investigating the measured shear for B- 
modes (Crittenden et al. 2002; hereafter C02), and most 
of the more recent cosmic shear studies perform such an 
analysis. Gravitational lensing effects cause a shear field 
which, to leading order, only contains E-modes. The fact 
that early cosmic shear measurements contained a signifi- 
cant B-mode contribution of the shear was interpreted as 
remaining systematics. As has been pointed out by several 
groups (e.g., Heavens et al. 2000, Crittenden et al. 2001; 
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Croft & Metzler 2000; Catelan et al. 2000; Jing 2002), 
B-modes in the observed shear field can be caused by in- 
trinsic alignments of source galaxies. This effect can be 
controlled in principle if sufficiently accurate photometric 
redshift information about the source galaxies are avail- 
able (King & Schneider 2002, 2003; Heymans & Heavens 
2003), but this is not the case for presently available sur- 
veys. In addition, B-modes are caused by gravitational 
lensing due to the redshift clustering of source galaxies, 
but the resulting signal is very small (Schneider et al. 2002; 
hereafter S02). If there are inhomogeneous systematic ef- 
fects, such as a spatially dependent shear calibration er- 
ror, B-modes would also be generated (Guzik & Bernstein 
2005). 

The separation of the shear field into E- and B-modes 
is usually performed at the level of second-order statis- 
tics. Two standard methods for this have been estab- 
lished in the literature. One of them is the aperture dis- 
persion statistics (Schneider et al. 1998) which is most 
conveniently expressed as an integral over the shear two- 
point correlation functions £± (see C02; S02). The second 
method is the calculation of shear correlation functions 
which contain only E- or B-modes (C02, S02). As was re- 
cently shown by Kilbinger et al. (2006), these two meth- 
ods suffer from the fundamental problem that in order 
to apply them, one would in principle need to measure 
the shear correlation functions to zero separation (in the 
case of the aperture dispersion) or to infinite separation 
(for the E/B-mode correlation). Both of these measure- 
ments are impossible, as pairs of galaxy images blend if 
they have too low an angular separation, and the finite 
size of observational fields prevents measurements at ar- 
bitrary separations. This problem of principle becomes a 
problem in practice if one studies the presence of B-modes 
on small and large angular scales. 

In this paper, we present methods for overcoming this 
problem of principle, by defining integrals over the shear 
correlation functions which (i) extend over the finite range 
"$min < i? < i^max where the range of integration can be 
chosen arbitrarily, and which (ii) provide a clean separa- 
tion of the shear signal into E- and B-modes. The first set 
of such integrals is termed the circle statistics; it is de- 
fined as the correlation of the mean tangential and cross- 
component of the shear on two circles with different ra- 
dius. As we will show in SectJH this correlation can be 
written in terms of integrals over the shear correlation 
function. Since the mean tangential (cross) shear on a cir- 
cle depends solely on the E- (B-) mode of the shear field, 
this provides a clean separation of both modes. In Sect. 01 
we consider the correlation of the mean tangential and 
cross shear over two rings which can as well be written in 
terms of integrals over the shear correlation functions £± . 
The relation of these new statistics to the aperture mass is 
briefly discussed in Sect.0J We then present a generaliza- 
tion of the circle and ring statistics in Sect.[S]which allows 
the choice of very general weighted integrals over the shear 
correlation functions that provide clear separation into E- 
and B-modes. We briefly conclude in Sect.EJ 



2. The circle statistics 

2.1. Definition 

Consider a circle of radius 9 around the origin and define 
the mean of the tangential and cross component of the 
shear on this circle, Ct and C x , respectively, measured with 
respect to the center, 

C(9) = C t (8)+iC x (6) 

= d<P (7t+i7x)(^) ■ (1) 

The tangential and cross-components of the shear are re- 
lated to the Cartesian shear components 7° = 7i + 172 
through a rotation, 

7t+i7x --7 c e" 2iv , (2) 

where ip describes the polar angle on the circle. As is well 
known (see, e.g., C02, S02), C t is sensitive only to the E- 
mode of the shear field, whereas C x is sensitive only to the 
B-mode. 

We next consider the correlation between the mean 
shear on two concentric circles with radii 9\ and 82, 

(C(*0C(fe)> = ^ r ^ e-*C*+"0 
Jo 2lT Jo 27r 

x { 1 c {e 1 ^ 1 ) 1 c {8 2 ^ 2 )) , (3) 

where we made use of P)l. Correspondingly, we define the 
correlator 

Jo 27r Jo 2lT 

x (7 c (0 1 , ( p 1 ) 7 c *(0 2 ,^ 2 )) , (4) 

where the asterisk denotes complex conjugation. 
Expanding the correlators then yields 

(C(0i)C(0 a )) = (Ct(0i)C t (9 2 )) - (C x (0!)C X (0 2 )) , (5) 

and the imaginary part vanishes if we assume that the 
shear field is parity invariant (Schneider 2003); this as- 
sumption is equivalent to the requirement that the mixed 
shear correlator (7t7x) = 0. Similarly, 

(C(9 1 )C*(6 2 )) = (C t (^)C t (e 2 )) + (C X (9 1 )C X (6 2 )) , (6) 

so that the separation into E- and B-modcs is achieved 
with 

(Ct(^)Ct(fe)) = \ [{C(6 1 )C*(9 2 )) + (C(0!)C(0 2 ))] , (7) 
<C X {6 X )C X (9 2 )) = ~ [(C(6 1 )C*(9 2 )) (C(8 1 )C(e 2 ))] . (8) 

2.2. Relation to the shear correlation functions 

The shear correlators that occur in © and J3J can be 
expressed in terms of the usual shear correlation function 
£± by rotating the Cartesian shear into tangential and 
cross components relative to the separation direction of 
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the two points with polar coordinates (6%, tpi) and (02, <p 2 ). 
We then hnd that 



< 7 c (0i,^i)7 c (02,^)) = e 4l ^-W , 

( 7 c (^^i)7 c *(^^2)) = e+w, 



(9) 
(10) 



where $ and ip are the polar coordinates of the separation 
vector, 



$e lv = 2 c i( 



More explicitly, the separation is 

d 2 = e\ + el - 26162 cos(A^) , 



(ii) 



(12) 



where Alp = tp2 — ipi is the difference of the polar angle 
of the two points under consideration. In the calculation 
of 121 we need the square of the expression 



e 2i v e -i(vi+ i P2) 



02 e -iAv, + 2 e iA v _ 2 e x 6 2 
T 2 



(13) 



The integrand in @ then depends only on the difference 
Aip of the polar angles, so that one of the angular integrals 
can be carried out. Furthermore, we note that the integral 
is real, since the imaginary part of the integrand is an odd 
function of Aip. This then allows us to write the correlator 
of the circle statistics as 

r7I dAp £_(tf) 



<C(0i)C(0 2 )) = 



o 4 



[(0 4 + 6 A 2 ) cos(2A^) 



-401 02(0? + 2 ) cos(A<^>) 



(14) 



We can now transform the integral into one over the sep- 
aration i9. Using 112[) we hnd for the transformation of 
variables 



dA<p = 



2r) dfi 



VWi + e 2 y - tf 2 ^ 2 - (#2 - e 1 y 



(15) 



The ^-integration then extends from 9 2 — B\ to 0i +02, 
where we set 0i < 02 without loss of generality. Then, the 
final expression can be written as 



(C(0i)C(0 2 )) 



102-9 

where the dimensionless function Yl is 
Y-(x,r,) ' 



(16) 



irn 2 x 2 a/(1 + rj) 2 — x 2 ^J x 2 — (1 — rf) 2 
x [(1 - x 2 ) 2 - 2 V 2 (2 - x 2 ) + ?7 4 (6 + 2a; 2 + x 4 ) 
- 2 V 6 (2 + x 2 )+ V s ] . (17) 
In a similar way, we obtain for the other correlator Q 

<C(0i)C*(0 2 )) = e+ W cos(2A^) 



' 1+tf2 dv) / § 61 , 



?2 <?2 



with 



Y+{x,rf) 



x 2 [(1 ~ x 2 ) 2 - 2r] 2 x 2 + f] 4 ] 
r\ 2 i^\J (1 + rf) 2 — x 2 \J x 2 — (1 — rf) 2 



(19) 



We therefore get for the E-mode of the circle statistics 



(Ct(0i)C t (0 2 )> 



10-2-8! 



di) 
2!) 



i) 



0i 

02 ' 02 

and the corresponding B-mode reads 



(20) 



(C x (0i)C x (0 2 )> 



1+02 



d$ 
2d 



(21) 



Hence, we have obtained a statistics which is able to 
perform the E-/B-mode separation of the shear correla- 
tion functions on a finite interval in angular separation, 
02 — 0i < i9 < 0i + 02- In particular, it avoids the use of 
correlation functions towards zero separation which, as we 
argued before, are difficult or impossible to be obtained 
from observations and whose neglect causes artificial E- 
/B-mode mixing in the aperture statistics. 

However, the circle statistics is not very practical, since 
the functions Y± have a i? -1 / 2 singularity at both integra- 
tion boundaries. This singularity has a purely geometrical 
origin, in that at the minimum and maximum separation 
of points on the two circles, the circles are 'mutually par- 
allel', so that pairs with these two extremum separations 
are particularly frequent. Applying the foregoing equa- 
tions to observed correlation functions with their noise 
will therefore yield estimates for the circle statistics which 
will be dominated by the noise just at the boundaries. It 
it is therefore desirable to find an expression for the E- 
/B-mode separation on a finite interval which does not 
contain such a singularity. This will be provided by the 
ring statistics, which is discussed in Sect. OH 

2.3. Relation to the power spectrum 

If we assume that the shear field is caused by gravitational 
lensing, no B-modes are present. The shear field then is 
a pure E-mode field which derives from the underlying 
surface mass density k (we use standard lensing notation 
here). The power spectrum P(£) of the k field then speci- 
fies the second-order statistics of the shear field, in partic- 
ular the correlation functions. We have (see, e.g., Kaiser 
1992) 

e±W= r^P{t)\A{&i), (22) 

where the J„ are Bessel functions of the first kind. 
Inserting this expression into the E-mode circle statistics 
(H3J yields 

dil 
2tt 

where 



(C t (0i)C t (0 2 )) = 



P(l)V E I62, 



(23) 



• 1+ *> dx 



^ [J (vx)Y+(x,r]) + J 4 (z/x)Yl(x,77)] 
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Fig. 1. The filter function Ve relating the power spectrum 
to the circle statistics, see Eq. <|23|) . for three different 
choices of the ratio of radii r/ — 9\j9 2 . 

Hence, the function Ve(^, rj) describes the filter function 
with which the circle statistics is related to the lensing 
power spectrum; it depends on the ratio r\ of the two radii. 
The filter function is plotted in Fig.^for three values of 
rj. It shows an oscillating behavior which reflects the os- 
cillations of the Bessel functions. 

The behavior of Ve for small and large values of v 
can be obtained analytically. Considering small v first, we 
can expand the Bessel functions in the foregoing equation. 
The leading-order term in J4 is (yx) 4 , whereas the expan- 
sion of Jo contains terms proportional to (yx)° and (yx) 2 . 
However, we find that 



r 1+v dr f 1+rt c\r 

so that the leading-order term from the Jo part is also 
cx . Together, we conclude that Ve(^, rj) °c v for small 
is, as can also be seen in Fig.^ 

In the opposite case of large v, the behavior of the 
integral for Ve(i / , rj) is dominated by the 1/ y^-singularity 
of the functions Y± at both integration limits. Hence, the 
asymptotic behavior of Ve(' / , rj) can be obtained from an 
integral with very similar properties, 



da; 



Jq(^e) 

V x q ~ x 



and 



da; 



iiiyx) 



V4 



2 / VX 



^12 ( VX °\ 



Thus, for large v we expect Ve(^, rj) to be an oscillating 
function, where the amplitude of the oscillation decreases 
as v . Again, this behavior is seen in Fig.^ 
The vanishing of the B-mode implies that 

" 1+n dx 

'vx)Y + {x,rj) - J 4 (^a;)y_(a;,77)] 



f 1+ri dx 
VbM = / r- [Jo( 



Fig. 2. The functions Z + (thick curves) and Z_ (thin 
curves) as functions of the angular separation scaled by 
the maximum separation, x — f?/i9 m ax, f° r three different 
values of the ratio r\ — ^min/^max- The specific choice for 
the ring boundaries (|28J) and weight functions ijHS} was 
assumed. 

for all values of v and arbitrary < 77 < 1. As a conse- 
quence, in the absence of B-modes, the circle statistics can 
be calculated from £ + and £_ separately, since both terms 
occurring in (|2U[I are then equal. 

3. The ring statistics 

Instead of the shear on a circle, we consider here the shear 
inside a ring and define 



n = Kt+\Tl, 



C2 



d9 W(9)C{9) 



(24) 



where the integral extends over the annulus £1 < 9 < 
£2, and W(9) is a weight function which we choose to be 
normalized to unity, 



C2 



d9 W{6) = 1 



(25) 



£1 



As was true for the circle statistics, the real and imagi- 
nary parts of 1Z are solely due to E-modes and B-modes, 
respectively. We now consider the correlator of two such 
concentric rings, one with (1 < 9 < ( 2 as before, the other 
with £3 < 9 < £4. If d < Q for all i < j, as we shall assume 
from here on, the rings are non-overlapping. The maxi- 
mum and minimum separation of points from one ring to 
the other are $ max = C2+C4 = * and tf min = C3-C2 = 
respectively. The correlator reads 



(KK) 







C2 d9 1 Wi(6i) f* d9 2 W 2 (9 2 ) {C{9 1 )C{9 2 )) 
d(9i Wx{9i) J d9 2 W 2 (9 2 ) (26) 
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Fig. 3. The filter function We relating the power spectrum to the ring statistics, see eq. <|32|) . for three different choices 
of the ratio of radii rj = 61/62. Dotted curves indicate negative values. 



Interchanging the order of integration, this correlator can 
be written in the form 



f 1 dx 

{ nn) = / — z_ (a, »?) 



(27) 



where 77 = i9 m i n /i9 m ax is the ratio of minimum and maxi- 
mum separation of points between the two rings, and the 
function Z_(x, 77) can be calculated explicitly by changing 
the integration order in 126|) . In general, the resulting form 
of the function Z_ will be complicated, due to the numer- 
ous case distinctions in the range of integrations. The task 
simplifies considerably if special ratios of the ring radii are 
chosen. In the appendix, we provide explicit formulae for 
Z_ for the following choice: 



Ci 
Cs 



1 _ V -r A 3(1 — 77) T 

* ; C 2 = n - V 



8 

5/7 + 3 
8 



T /• 3?7 + 5 T 
* ; C4 = -4^* • 



(28) 



This choice still allows the freedom to choose the range 
77 \& < ?9 < ^ over which the correlation function is probed 
in the ring statistics. In analogy, the other correlator in- 
volving reads 



(nn*) 



dx 



Z + (xV)Z + (x,r,) 



(29) 



To be specific, we will choose for the weight functions 

30(^-Ci) 2 (C 2 -ei) 2 



Wi 



(C2 - C1) 5 



W 2 (6 2 ) 



3O(0 2 -C 3 ) 2 (C4 



(30) 



(& - Cs) 5 

In Fig. Elwe have plotted the functions Z± for three values 
of 77. Both functions and Z- have two roots in the 
interval rj < x < 1, starting off with positive values near 
the edges of their range of support. We will explain this 
qualitative behavior in Sect. 15.21 below. 

The E- and B-modes of the ring statistics are given by 
the combinations 

<^) E = \ «ftrc*> + (mi)) , 

= J —{U(^)Z+(x, V )+^(x^)Z_(x,r 1 )], 

(31) 

^(x^)Z_(x, V )\ . 



(KK) 3 = 



~ {(nn*) - (Tin)) 



In case the shear field derives from gravitational lensing, 
its B-modes vanish, and the ring statistics (R.Tl) E can be 
expressed in terms of the power spectrum P(£) of the pro- 
jected surface mass density. In analogy to the circle statis- 
tics, we obtain 



Jo 2vr 

■re 

f 1 dx 

(^V)= J 1^ \- J o(vx)Z+(x,r]) + J A (vx)Z-(x,r))] 



(32) 



where 
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Fig. 4. The ring statistics (TZTZ) E (|32l) as a function of the 
maximum separation W = $ max for which the correlation 
functions need to be known, for three different ratios r\ = 

max- 



In Fig. we have plotted the function We{^,v) which is 
the filter function that relates the power spectrum to the 
ring statistics, for three different values of 77. As we can see, 
We is a smooth function for small v, whereas it strongly 
oscillates for larger values of v. Asymptotically, We(^, 77) 
goes as v A for small v, and the amplitude of oscillations 
falls off as v~ 7 for large v. 

The ring statistics (1Z1Z) E is plotted in Fig. for a 
ACDM model with erg = 0.85 and a mean redshift of 
source galaxies of z w 1.5. For this calculation, the non- 
linear power spectrum according to the prescription of 
Smith et al. (2003) has been used. We see that (1Z1Z) E 
has a very small amplitude which decreases for larger 77. 
Hence, the narrower the interval is over which the correla- 
tion functions are considered, the smaller is the measured 
shear signal. For all values of 77, the signal of (1Z1Z) E is 
much smaller than the characteristic values of the shear 
correlation functions. The mathematical reason for this 
can be seen in Fig.[5] in combination with (|31|1 : to derive 
(1Z1Z) E , the shear correlation is integrated over the func- 
tions Z± which has two roots in their range of support. 
Another way to see this is by considering the function 
We (Fig.EJ in connection with (|32(1 . Due to its oscillat- 
ing behavior, the resulting integral in (|32[l is expected to 
be small. The narrower the ring is (i.e., the larger 77), the 
more oscillations does the function We show near its broad 
maximum, and the larger is the near-cancellation of power 
in the integral l|32|l . 

The vanishing of the B-mode then implies that 

f 1 dx 

Wb(v,v) = J — [3 (vx)Z + (x,r]) - J 4 (vx)Z-(x 1 i])} 







(33) 



for all values of v. We have checked numerically that the 
functions Z± indeed satisfy this relation. 



4. Relation to the aperture statistics 

The definition (24) of the ring statistics is very similar to 
that of the aperture mass (Kaiser 1995; Schneider 1996), 
which is defined as 



M, 



d 2 7?K(tf)[/(|tf|) , 



(34) 



where U is a radial weight function and the integral ex- 
tends over the support of U. Provided U is a compensated 
filter function, / di?0i/(i?) = 
terms of the tangential shear, 



M ap can be written in 



M, 



d 2 7? 7t (tf)Q(|tf|) 



where the weight function Q is related to U by 



(35) 



(36) 



The definition of 7Z is equivalent to l|35(l in the absence of 
B-modes, if we set W(8) = 2ir8Q(9). It is therefore inter- 
esting to see whether the ring statistics can be expressed 
also in terms of the surface mass density. This indeed is the 
case, since l|5rj|) can be reversed. Multiplying (I3rj)l with 9 2 
and differentiating with respect to 9 yields the differential 
equation 



U' 



-Q'-- n Q, 



(37) 



which can be readily integrated, using the boundary con- 
dition U{9) — >• as -*• 00: 



U{6) 



2dtf 



Q(tf) - Q{8) 



(38) 



It is straightforward to check that U is a compensated 
filter function. Furthermore, if Q(9) = for 9 > 6* max , 
then also U(9) — for 9 > 9 max . Therefore, there is an 
equivalent [/-filter with which the ring statistics can be 
expressed in terms of the surface mass density; U is a posi- 
tive constant for < 9 < d , then decreases and eventually 
becomes negative somewhere in the interval £1 < 9 < (2, 
and gets zero again for 9 > £2- 

5. Generalization 

The relations (|31|l that define integrals over the shear cor- 
relation functions which are sensitive separately to E- and 
B-modcs are very similar to analogous equations for the 
aperture dispersions (see S02, C02). We could easily find 
other forms of the function Z± by taking different choices 
of the weight functions Wj. 2 entering the definition of the 
ring statistics. 

Therefore, the question arises whether one can find a 
wide class of functions, analogous to Z± in (|31|l which 
provides a separation of the shear correlation into E- and 
B-modes. This question will be answered here. 
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5.1. The general E/B-mode decomposition 

Thus, we define second-order shear statistics in the follow- 
ing form: 



E = 



D 



dtfi?K + (tf)T + (tf)+£_(tf)T_(tf)] , 
dtftf[e+(tf)T+(tf)-£_(tf)T_(tf)] , 



and look for pairs of function T± ($) for which E contains 
only E-modes, and B contains only B-modes. In terms 
of the E/B-mode power spectra, the correlation functions 
read (see S02) 



die 
dee 



Jo [PbW+PbV)] 



Jo ^ 



(40) 



Inserting (|4(j|) into 1|39|) and interchanging the order of 
integration, we obtain 



E 



dee 

2tt 



p E (e)[t + (e) + ue)} 



+p B (£)[t + (£)-t.(£)}) , 
r°° de e / 

B = J — {P E (£)[t + (£)-t.(£)} 

+PB(e)[t+(e)+t-(e)}) , 



(41) 



where we have defined the Hankel transform of the func- 
tions T±, 



t+(e) = / dt?<?r + (0)Jo(&?) 

Jo 



Ue) 



dtftfT_(tf) J 4 (&?) 



(42) 



Thus, the functions E and B contain only E- and B-modes, 
respectively, if and only if the two Hankel transforms are 
identical, 



t+(£) = t-(£). 



(43) 



Any pair of functions T± which satisfies l|43[) leads to a 
separation into E- and B-modes. The functions Y± for the 
circle statistics and the functions Z± for the ring statistics 
are special examples of this, as are the pairs of functions 
T± for the aperture statistics and S± for the shear disper- 
sion, as defined in S02. 

We will show next, by explicit construction, that one 
can choose one of the two functions, T + or T_, arbitrarily 
and the other function can be uniquely calculated such 
that (|43[) is satisfied. First we assume that the function T_ 
is chosen; then we calculate T + from the inverse Hankel 
transform, 



T+(0) = / dt£J Q (£0)t + (0) . 
Jo 



Due to the requirement l|43|) we can replace t + by t- and 
obtain 



<M£J {M) / d0 0J 4 (*0)T_(0) 



d6 eT_{9)G{$,9) , 



(39) 

where we defined the function 



G{&,0) = / d££J (M) J 4 
Jo 

In S02 it was shown that 



(44) 



(45) 



f 4 1 2il 2 \ 1 
G (^> °) = {02 ~ —J - *) + fni ~ *) > (46) 

where H(x) is the Heaviside step function and <5d the Dirac 
delta function. Thus, 



T+(tf) - T_(0) + d6 6T_{6) (± 



4 12tf 5 



4 



(47) 



Hence, for any function T_(i?) we can calculate from (|47|l 
the corresponding function T + so that this pair of func- 
tions satisfies <|43[1 . In a similar way, T_ can be obtained 
from a chosen function T + , from 



r_(tf) = r+(t?) + jT d0 0T + (0) (^- 



4 120 2 



(48) 



5.2. General E/B-mode decomposition on a finite 
interval 

We now return to the original question of this paper, 
namely to find a separation of E- and B-modes from corre- 
lation functions given on a finite interval. Thus we choose 
a function T_(i9) which is non-zero only on the interval 
$min < # < $max- The upper limit of the integral in l|4T|l 
can then be replaced by i3 max , and we see immediately 
that T+(i?) = for i9 > tf max . If we require that T+(tf) 
vanishes for i? < $ m in, then the following two conditions 
need to be satisfied: 



(49) 



Thus, if we choose a function TL with is non-zero only on 
a finite interval, the corresponding function T + vanishes 
outside this interval as well, provided T_ satisfied the two 
integral constraints (|49|l . The argument which led to the 
constraints (|4l?|> only applies for # min > 0; for $ min = 0, 
these conditions become obsolete. 

Conversely, if we choose a function T + which vanishes 
outside the interval $ m i n < ^ < $max) the corresponding 
function T_ is zero outside the same interval if T+ satisfies 
the two constraints 



d$-dT + (d) = 



di?tf 3 T+(tf) 



(50) 



if tf E 



< 00. 
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In fact, it is straightforward to show that the defini- 
tion of T + in terms of T_ as given by 147|) implies the 
relations l|5(J|l provided i? m i n > 0; this latter condition is 
needed to guarantee the existence of the integrals that oc- 
cur. Similarly, the definition (|48|l implies the validity of 
(113) if tfmax < oo. 

We can check these results in a different way, by ex- 
panding the Bessel functions into a power series and com- 
paring the result, when inserted into H42|) term- by-term. 
Using 



Jo(z) 



z 1 



(-z 2 /4) A 



^ o (fc + 2)!(fc + 2)! 



\ " 



fe=0 v 



(-z 2 /4) fc 
"4)1 



(51) 



we see that the integral (|42() of each term is guaranteed 
to exist if T± have finite support, which we assume hence- 
forth. The first two terms of Jo, when inserted into the 
first of (|42|) . vanish because of ||SDJ. We next compare 
terms proportional to ^( 2fc + 4 ) of t±. For i_, this term is 



d-d dT-(tf) 



(_^2 /4) (fe+2) 



fc!(fc + 4)! ' 
whereas for t + , the corresponding term reads 

t +k = dWT + ({))- [ ' 



(fc + 2)!(fc + 2)! ' 
Inserting now the expression (|47|l for T + , we find 

(_4)-(fc+2) 



(52) 



(53) 



t+k 



(fc + 2)!(fc + 2)! J Q 
T_(i?)+ / d9 8T_(0) 



dtf tf( 2fe + 5 ) 



4 
tf 2 



12# 



2 \ 1 



(54) 



For the second term, we interchange the order of inte- 
gration, after which the ^-integral can be evaluated. This 
then leads to 



x T_(i?) 



1 



(_^2 /4) (fc+2) 

(fc + 2)!(fc + 2)! 
2 6 



fc + 3 



(55) 



Using the definition of the faculty, it is now easy to show 
that t + k = t~k for all k. Hence, for weight functions with 
finite support, the equivalence (|43(l can be demonstrated 
for each order of I. 

It is interesting to consider two special cases of the 
definitions E and B that have been investigated in the 
literature. The first case is that of the aperture dispersion, 
as discussed in S02. The corresponding functions there 
were also denoted by T±, 



T + {x) 

x\J A — x 
+ IOOtt 



6(2- 15a; 2 ) 



1 



HI) 



(120 + 2320a; 2 - 754a; 4 + 132x 6 - 9a; 8 ) (56) 



for < x < 2, and zero otherwise, and 



T_(x) 



192 
35^' 



x' 

T 



2\ 7/2 



H(2 - x) 



(57) 



Here, x = $ /9q, where 9q denotes the aperture radius. The 
functions are plotted in S02; they both vanish for x > 2. 
If we consider T + as given, then T_ follows from (|48|) : 
since T + satisfies l|50fl. T_ shares the property of T + to 
vanish for x > 2. Conversely, consider T_ as given. Since 
T_ = for x > 2 one obtains immediately from H47JI that 
T + shares this property. However, T_ does not satisfy the 
relations (|49|l : in fact, since T_ oc x for x <^1, the second 
of the integrals in (|49|l does not exist. 

The second example is that corresponding to the shear 
dispersion in a circle, where the corresponding weight 
functions have been denoted as S± in S02. They are 



S+(x) 



l r 



4 arccos 



(?W 



4-rr 2 



H(2 



(58) 



_ . , a;V4^ 2 (6-a; 2 )-8(3-a; 2 )arcsin(x/2) 
S-(x) = 59 

for x < 2, and S-(x) = 4(x 2 - i)/x i for x > 2. Thus, 
whereas S+ is confined to a finite interval < x < 2, S- 
is not. The reason for this is that 5+ does not satisfy (ISOJ); 
in fact, S+(x) is a non-negative function. 

The constraints Q49JI and H5()|l also explain the quali- 
tative behavior of the functions Z± shown in Fig.0 The 
two constraints in either case require that the functions Z± 
have two roots in their range of support. As this was the 
reason for the small amplitude of the resulting shear signal 
shown in Fig.0] it becomes clear that this small amplitude 
is not due to the specific choice of our ring statistics, but 
is a generic feature of the more general statistic E. 

6. Summary and Conclusions 

In this paper we have constructed new second-order cos- 
mic shear statistics which can be used to separate E- and 
B-modes from measurements of the shear correlation func- 
tion £± on a finite interval. The circle and ring statis- 
tics were constructed, based on the fact that the cross- 
component of the shear averaged over a circle vanishes for 
pure E-mode shear, and the tangential shear component, 
when integrated over a circle, vanishes for a pure B-mode 
shear. Both of these statistics are obtained by integrating 
the shear correlation functions, multiplied with properly 
defined weight functions, over a finite interval in angular 
separation. By choosing the parameters of the ring and/or 
circle, the interval on which the shear correlation shall be 
probed can be adjusted. 

Motivated by the result, we have shown that there ex- 
ists a very large family of such weight functions which can 
be used to separate E- and B-modes on the second-order 
statistics level from measured shear correlations on a finite 
interval. Hence, there is a large freedom to choose these 
weight functions, of which the circle and ring statistics are 
special cases. 
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The resulting signal of the E- and B-mode shear statis- 
tics is small; this was explicitly demonstrated for the ring 
statistics (see Fig.0J, but as we argued in Sect . 15. 21 such a 
small signal is expected in general and can be traced back 
to the integral constraints (|49|l and l|50[l that the weight 
function T± have to satisfy if they are to probe the shear 
correlation for E-/B- modes over a finite interval in angu- 
lar separation. The small signal implies a correspondingly 
small signal-to-noise for these statistics. In particular, if 
fairly narrow ranges of ^ m i n < $ < i? max are to be probed, 
the corresponding shear data must have high statistical 
accuracy. We therefore expect our results here to be of 
particular relevance for the upcoming large cosmic shear 
surveys. 

The ring and circle statistics presented here can be 
generalized to higher-order cosmic shear statistics as well. 
The aperture statistics that were defined to quantify third- 
order shear statistics (Jarvis et al. 2004, Schneider ct al. 
2005) are plagued with the same problem as the second- 
order aperture statistics pointed out by Kilbinger et al. 
(2006). Hence, one could generalize the circle statistics to 
(C(8i)C(92)C(9 3 )) and related expressions (by taking com- 
plex conjugates of different factors), project out the E- 
and B-modes, and express the result in terms of the shear 
three-point correlation functions. By taking weighted inte- 
grals over finite ranges of 6i , the third-order ring statistics 
can be obtained. In this way, also the third-order shear 
can be tested for the presence of B-modes. 
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Appendix A: Explicit calculation of the functions 

Z± 

In this appendix, we carry out the steps that lead to the 
form of 127|) , and provide explicit expressions for the func- 
tion Z_. Starting with the final form of l|26l) . then first 
interchanging the order of integration between 9 2 and "d, 
after that changing the order of integration between 9\ 
and i9, we obtain 



(UK) = Ve-W / d6 1 W 1 (9 1 ) 

x / d6 2 W 2 (e 2 )Y(-,-±) . (A.l) 
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<S> 



■C 



Cl 

Ci 



C2 = C3 = C 4 

C 2 = x + yi 



5t/ + 3 



C3 — C4 — C5 — Ci 



c 6 = a; - yi 

3?j + 5 



and the functions Wi(yi) — ^>Wi(yi^). For the choice l|30[l 
of the functions W, , we obtain 



^i(yi) 

W 2 (y 2 ) 



15(8 yi + 77 - 1) 2 (3 - 377 - St/Q- 

2(1-7^ 
15(8y 2 - 577 - 3) 2 (3t7 + 5 - 8t/ 2 ; 
2(1-77)5 



(A.4) 



The expressions for (1ZH*) are identical, except that £_ 
and Y_ are replaced by £ + and Y + , respectively. From 
these explicit expressions, the functions Z±(x,rj) can be 
read off readily. They are plotted in Fig. 



0i 



Fig. A.l. The trapezium (solid line) is the integration 
area in the 0i-7?-plane for Eq. (|A.1(1 . The devision into six 
triangular regions is indicated, the numbers correspond to 
the indices of the integration bounds in (|A.3(1 . 



The integration range in the 9\-"d plane is sketched in 
Fig. IA.1I For this figure, we have chosen the values of 
Ci to satisfy the relations (|28|) : as it becomes clear from 
this figure, the integral then naturally splits up into six 
parts. Changing to dimensionless variables by defining 

= xtf ; 01 = yi* ; 2 = , (A.2) 

with ^ = (2 + Ci> we can then write 

■ i J at x Jbi 



i=l 



bi 

x [ 'dy 2 w 2 (y 2 )Y_ (-A) , (A.3) 

where the various integration bounds are given as 

3t/ + 1 



a\ = r? ; A 6 = 1 ; A\ = a 2 = a 3 = 



4 



A A 1 + 7? a A ?/ + 3 

a 4 = a 5 = A 2 = A 3 = — — ; a 6 = A 4 = A 5 = — - — 



61 = -64 = -y3 5 = — x 

Bi = B 3 = B± = B§ = - - ^ 

b 2 =h = — 
B 2 = o 3 = -6 6 = — x 



